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ABSTRACT 



A control system for computing economic dispatch 
signals for each operating area within an intercon- 
nected power system, utilizing an algorithm for solving 
the interarea coordination equations resulting in the 
definition of a common reference running cost for the 
interconnection, and enabling the explicit solving of 
the individual operating area running costs, thus 
avoiding any need for an iterative solution. The com- 
putational burden is shared between a system com- 
puter which makes periodic power and/or running 
cost assignments for each operating area, and operat- 
ing area computers which independently calculate 
specific dispatch signals for the generators within their 
areas. 

10 Claims, 5 Drawing Figures 
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ECONOMIC DISPATCH TECHNIQUE FOR 
INTERCONNECTED POWER SYSTEMS 
This is a continuation of application Ser. No. 
293,010, filed Sept. 28, 1972, now abandoned. 

BACKGROUND OF THE INVENTION 

A. Field of the Invention 

This invention lies in the field of computer controlled 
interconnected power systems and more particularly in 
the field of interconnected power systems wherein 
there is provided explicit periodic calculation of dis- 
patch signals for each operating area within the inter- 
connected system, based upon the tie line flow of 
power between the areas and the power generation 
within each such area. 

B. Description of the Prior Art 

Electric power generating systems generally are com- 
prised of a plurality of generating units of differing effi- 
ciencies and having differing absolute and incremental 
costs of power generated. Generally speaking, systems 
are comprised of operating areas, each of which is a 
generating system network within the larger system. An 
interconnection, or pool, as used herein, is a group of 
such discrete operating areas, interconnected for eco- 
nomic operation. A control area, as used herein, is an 
operating area or interconnection of operating areas 
which maintains the net flow of power across its bound- 
aries at or near a scheduled value which is usually re- 
vised periodically. Depending upon the contractual re- 
lationship, an interconnection may comprise one or a 
number of constituent control areas. Within any given 
control area, generating units are regulated by a load 
dispatch system so as to match system generation with 
system load on an economic basis. 

In the operation of such interconnected systems, it is 
of course desirable to optimize system stability and 
minimize uneconomic response, particularly from non- 
regulatory units. A copending application filed by the 
same inventor, U.S. Ser. No. 163,894, filed July 19, 
1971 and titled "COMPUTER CONTROLLED CO- 
ORDINATION OF REGULATION AND ECONOMIC 
DISPATCH IN POWER SYSTEMS," sets forth an im- 
proved control system for achieving the above objec- 
tives. In both prior systems and the above improved 
system, frequency and tie line power flow for a control 
area are continuously monitored and compared with 
scheduled frequency and tie line flow, respective differ- 
ences being combined to produce an error signal, 
called "area control error" or ACE signal, representing 
the difference between generation and load. In the im- 
proved control system, the ACE signals are combined 
with a summation of telemetered actual generation sig- 
nals and supplied to an incremental cost computer 
which generates therefrom a calculated cost signal that 
is converted to a generation control signal and trans- 
mitted to respective economic generating units. 

The incremental cost computer of the referenced sys- 
tem, which describes the operation of an interconnec- 
tion like the Pennsylvania-New Jersey-Maryland inter- 
connection, periodically makes a detailed dispatch cal- 
culation. This is an iterative calculation which refer- 
ences the incremental cost of every generator within 
the interconnected system to a common reference 
buss, the iteration being continued until all incremental 
costs are equal with respect to the reference buss. This 
calculation determines the incremental cost, or lambda 
(X), of delivered power for the entire system and for 
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each company (or economic unit) in the intercon- 
nected system, from which ratios of the respective 
company lambdas to the system lambda are also calcu- 
lated. The system lambda is then modified by the re- 

5 spective ratios before being transmitted to the member 
companies. Algorithms for the iterative dispatch calcu- 
lations are well known in the art and are in general use 
in power interconnection control systems. 
In the above prior art systems, each sub-area com- 

10 pany, after receiving a X signal, recalculates its own 
economic dispatch signals, on the basis of its own sys- 
tem model, not taking into account the entire system. 
Consequently, there develop discrepancies between 
such sub-area calculations which impair the accuracy 

15 of the economic dispatch. Also, as a result of the 
lengthy nature of the complete area dispatch, it may be 
expected that the calculation of the lambda ratios will 
be made relatively infrequently (every 5 or 10 minutes) 
because of the penalty in computer usage for essen- 

20 tially redundant calculations. In the interval between 
such calculations, the accuracy of the lambda ratios de- 
creases, due to changing system conditions. Further, in 
recalculating the economic dispatch functions for each 
sub-area, the non-linearity of the local incremental cost 

25 curve aggravates any error in the system lambda. 

Due to the low incremental cost and large size of 
most generating units, the system incremental cost 
curves are characterized by flat sections covering large 
blocks of energy combined with much steeper sections 

30 toward the upper end of the curve. At the same time, 
most steam units respond sluggishly at best to an incre- 
mental cost signal. This combination of circumstances 
makes it all but impossible for the system to reset, or 
compensate the previously calculated lambda value, 

35 and the presumed system lambda can rapidly become 
so inaccurate as to be unstable. 

It is to be noted that without considering transmission 
costs, all lambdas of the constituent companies, or sub- 
areas within the interconnected system, are made equal 

40 in order to minimize overall cost. However, because of 
transmission costs, the different companies do in fact 
have different lambdas at which they deliver power. In 
view of such transmission costs, it becomes desirable 
that the incremental cost of power delivered to the 

45 load, for any two or more generating units within the 
system, have the same incremental costs. However, for 
interconnected systems, the expense of power delivery 
to a boundary between sub-areas, e.g., from one com- 
. pany to another, represents the cost paid by the receiv- 

50 ing company, or the cost of delivering to an economic 
load. Thus, the controlling criterion should be that the 
interchange incremental power costs, being the deliv- 
ery costs, be the same for all sub-areas, and be mini- 
mized. This criterion also has the advantage of simplify- 
ing accounting procedures, since billing is done on the 
basis of average interchange power between compa- 
nies. The system and method as disclosed in this appli- 
cation adopt this criterion, and utilize a unique algo- 
rithm for explicit computer calculation of optimum sys- 
tem and sub-area running costs. 

SUMMARY OF THE INVENTION 

It is the primary object of this invention to provide a 
65 method and means for enabling computer controlled 
economic dispatch for an interconnected power system 
wherein the dispatch signals are solved for explicitly in 
terms of tie line flow between the sub-areas and power 
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generation within each area, and individual generation cost for the interconnection, allowing the explicit solu- 

signals for the separate control sub-areas (companies) tion of individual sub-area running costs. The proce- 

within the interconnection are calculated and transmit- dure for making the economic dispatch calculations is 

ted to such companies. as follows: 

It is another object of this invention to provide as 1. Each member company transmits, to a central 

computer controlled technique for, calculating eco- computer which makes the interconnection dispatch 

nomic dispatch for an interconnected power system calculation, at least that portion of its effective incre- 

which is quicker than prior art systems, optimizes com- mental cost curve within a specified band of its current 

puter usage, and enables more frequent calculation of load, in effect providing the interconnection with a lin- 

the dispatch so as to avoid the errors of slower prior art 10 earization of each company's effective incremental 

systems. , f cost function about such current operating point. 

It is another object of this invention to provide a sys- 2. The interconnection system computer adjusts each 
tern and method for minimizing the overall computa- member company cost curve by means of a compensa- 
tional burden in calculating economic dispatch for an tion factor for such company, which factor incorpo- 
interconnected power system, and in which the compu- 15 rates transmission cost data, and then combines the ad- 
tational burden is shared more evenly and without du- justed cost curves to get a total interconnection system 
plication between an interconnection computer and cost curve. The compensation factor for each company 
each of a plurality of company, or sub-area computers, is obtained by the unique algorithm which allows ex- 
with direct control of generation retained at the control plicit determination. The terms particular to each corn- 
area level. 20 pany which are required to calculate ithe compensation 
It is another object of this invention to provide a factors are transmitted from each such company to the 
computer controlled technique for determining eco- interconnection computer along with the incremental 
nomic dispatch for an interconnected power system cost curve information, at each periodic calculation of 
wherein the effect of variations in the slope of the ef- the economic dispatch signals, 
fective interconnected system incremental cost curve 25 .3, The interconnection system computer determines 
are minimized. an overall (reference) incremental cost, or system run- 
It is a further object of this invention to provide a ning cost, by comparing the current total interconnec- 
computer controlled technique for economic dispatch tion system load (as determined from the ACE and 
of an interconnected power system wherein inter- total generation signals) with the combined adjusted 
change incremental power costs are equalized and av- 30 system curve. The power (megawatt) load that should 
erage incremental interchange power cost is mini- be carried by each company is determined by compar- 
mized. ing the system running cost with the adjusted cost curve 
It is a further object of this invention to provide a for each such company, thus obtaining the assigned 
computer controlled economic dispatch technique for load. The assigned load is then transmitted to the com- 
an interconnected power system wherein overall pro- 35 pany, which then independently calculates the amount 
duction costs, including transmission costs, are mini- of load to be carried by each of its generators, 
mized, and which permits each member company to 4. For any company within the interconnection sys- 
follow its own dispatch procedures independently of tern whose dispatch system is predicted on a running 
the control system procedure. cost (lambda) signal, and which cannot accept a power 
It is a yet further object of this invention to provide 40 signal, the interconnection system computer deter- 
an economic dispatch technique for an interconnected mines the company lambda by dividing the system fun- 
power system wherein the average cost of power deliv- ning cost by the company compensation factor, and 
ered to inter-company (sub-area) ties is minimized, and transmits the desired lambda signal with the assurance 
such average cost is made available for accounting that the signal thus supplied reflects true running cost. 



charges between respective companies within the inter- 
connection. 



45 



BRIEF DESCRIPTION OF THE DRAWINGS 



It is still a futher object of this invention to provide FIG. 1 is an illustration of an interconnected system 
a technique for computing economic dispatch for an constituting one control area, and having three inter- 
interconnected power system, wherein there is ob- connected companies (sub-areas), with tie lines indi- 
tained a common reference running cost for the tie cated. 

connections of the system. • FIG. 2 is a block diagram which shows the flow of in- 

In accordance with the above objectives, there is pro- formation to and from the computers utilized in the sys- 

vided a computer controlled system of computing the tern of this invention, for an interconnection system 
economic dispatch for an interconnected power system 55 having three member companies, 

comprised of a plurality of sub-areas, and utilizing a so- FIG. 3 is a flow diagram illustrating the steps in calcu- 

lution algorithm for such computation which is in lating power assignments and running costs for compa- 

closed form rather than iterative form, thereby provid- nies within the interconnection system, 

ing an explicit solution of the dispatch algorithm based FIG. 4 is a flow diagram illustrating detailed steps in 
on recently observed values of tie line flow between 6Q calculating the participation factors and compensation 

sub-areas and power generation within each sub-area. factors. 

The novel algorithm of this system adopts the classical FIG. 5 is a block diagram illustrating the control of 
technique for minimizing a function (overall produc- generator units at the company level, 
tion cost) subject to specified constraints, except that 
additional constraint equations for the power flow in 
the tie lines between sub-areas are included in the deri- 
vation of the dispatch values. Such additional equations The discussion which follows will be based on the il- 
permit the definition of a common reference running lustration of a three-company pool, as shown in FIG. ,1. 



DESCRIPTION OF THE PREFERRED 
65 EMBODIMENT 
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The three companies are illustrated as companies A, B , c. Cost curve data. For the current generation of 

and C, there being two tie lines between company A each company, each computer 60 transmits a portion 

and company B, and a single tie line between company of its stored company cost curve centered about such 

A and company C. The first tie line between company current generation level. 

A and company B is from terminal a\ to terminal M 5 System computer 50 is illustrated to receive C.F. 

with the power flowing across such tie line having an data, generation data and cost curve data from the 

incremental cost lambda of Xj. The second tie line from companies, which data is receives periodically at a rate 

company A to company B is from terminal a2 to termi- which is optimized in view of system parameters, e.g., 

nal 62, having a lambda designated as X 2 . The tie line fate of change of load. The system computer makes pe- 

between company A and company C is between termi-. 10 riodic calculations of the running cost, or and com- 

nals a3 and cl, having a lambda designated as X 3 . pany running costs, and transmits either a running cost 

Several definitions are helpful at this point. Where signal or a corresponding megawatt (MW) signal, or 

the term computer is used, it is to be understood as both, to each company computer 60. Computer 60 in 

meaning a general purpose digital computer having a turn makes an independent calculation, based upon the 

stored program. Also, it is considered that standard in- 15 signal it has received, of dispatch signals which are as- 

out peripheral equipment is utilized, for receiving and signed to each generating unit withinits company. Such 

transmitting data to and from the computers. Also, as independent company calculations may be made in any 

used herein, "interconnected system" and "pool" are known way, and the specifics of such company calcula- 

synonymous, as are "company," "operating area," and tions as such do not form a part of this invention. It is 

"sub-area." Thus, an interconnected system, or pool, 20 important to note, however, that each company may 

refers to a plurality of companies, operating areas or use its own technique for dispatching power generation 

sub-areas which are interconnected by tie lines over assignments within its own boundaries, 

which power may be passed from one company to an- The classical solution to the basic (unitary area) dis- 

other. patch problem- is obtained by formulating economic 

Referring now to FIG. 2, there is shown a combined 25 dispatch as a< parameter optimization problem with 

block diagram and information flow diagram for the in- equality constraints. In the system of this invention, I 

terconnection system of this invention. A system com- take this same approach to the compound area (multi- 

puter 50, typically a general purpose digital computer company) problem, introducing only slight modifica- 

with a stored program and with facility for periodically tions into the usual mathematical developmenti but ar- 

receiving input data, carries out the mathematical func- 30 riving at a new algorithm. 

tions directed by the unique algorithm of this invention. In the course of the development hereinbelow, the 
Associated with each company there is a computer 60, following assumptions are made: : 
such that for the three company illustration there is a A. AH the functions used in the development are con- 
computer A associated with company A, a computer B tinuous functions of time. 

associated with company B, and a computer C associ- 35 B. In addition, the composite, effective incremental 
ated with company C. Each company computer has a cost functions of the (three) individual member corn- 
stored program for making independent dispatch cal- panies are monotonically increasing functions of total 
culations for controlling the generators within its com- generation: : 
pany, such calculation being based upon a megawatt or C; The only sources of power injection into any one 
running cost (X) signal calculated for it by the system 40 company area are that company's own generating units 
computer 50. Cost curve information for a given com- and the tie lines connecting it to other areas, 
pany is stored in its respective company computer, for D. The solution to the problem of system control 
periodic transmission to the system computer 50i It is (i.e., overall economic dispatch of the entire com- 
to be noted that such cost curve information must be pound area) comprises the solutions to the problems of 
periodically updated, as by re-calculation, to account 45 the several local control areas (i.e., the economic dis- 
for varying conditions (as when a generator is taken on patches of the several member companies), 
or. off line). Alternately, the cost curve data for each ' The solution is begun in the usual manner by defining 
member company may be stored and/or calculated in a function 
system computer 50, if it has sufficient capacity. The A x , f 
method of making such cost curve calculations is well H + r~ 
known in the art. where 

Each company computer is connected, through con- p represeiits ^ tota i cost 0 f generation 

ventional transmission links to specific data points of f represents a set of constraint equations and 

such company, as shown in more detail in FIG. 4, and 55 x represents a set of undetermined coefficients (the 

obtains periodically sample values of current load and - ^ multipliers), 

tie flows. Each company computer 60 periodically For ea ch area there is a constraint equation 
transmits to system computer 50 the following data: 

a. Compensation factor (C.F.) data, comprising tie /,-«* + «,+ ' ( s«mi-#©,-o 
flow information from the company, as well as the in- fiQ 

cremental transmission losses for the tie lines. Incre- where 

mental transmission loss data is stored in each company PD t represents the power delivered to the load in 

computer, and the incremental transmittal loss is deter- company /, 

mined for the then sampled value of each tie line, and PL, represents the transmission losses within corn- 
transmitted to the system computer. 65 pany /, 

b. Generation date. Each computer 60 makes a sum- PE } represents the power exported from company / 
mation of total company generation, which summation over one of the / tie lines having one terminal in 
is transmitted to system computer 50. company /, and 
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PGi represents the power generated within company 



(c) For Company A, 



= (SPG*}, 5 

for the m generating units in company 7. 

In addition, there are n additional constraint equa- 
tions for the n inter-member-company tie lines, 

/j-TO+f^^oj # k 
which constrains the power into one terminal of a tie- 
line to equal the power out of the' other .terminal. This 
set of constraint equations enables treatment of the 2n 
tie-line terminals as independent power sinks in the 
several companies. (In forming H it is convenient to 
add the /i and subtract the /j.) ■• 

With these constraint equations, the vector A be- 
comes, for the example case, 
V — (X,,,: V* Xjm ^jtt* X^j) ■ - 

Giyen assumption B, the conditions for economic dis- 
patch (which are the necessary conditions for a station- 
ary value of F subject to the indicated constraints) are 
provided by . ... 

A. The (n + 3 ) constraint equations themselves; 

B. The (three) partial derivatives of H with respect 
to PG«, PG fc , PG C . Of course, H could be written, in 
terms of, and differentiated with respect to, the 
power generated by all the individual machines 
within the interconnection, thus yielding simulta- 
neously the conventional conditions for economic 
dispatch within each company; and 

C. The (2/i) partial derivatives of H with respect to 
the power delivered to. the several tie terminals. 

There are thus (n + 3) + 3 + 2n ~ 3« 4- 6. equations 
with 3n + 6 unknowns. It is important to note, since tie 
flows are treated as independent variables, in differen- 
tiating H with respect to generation, that under as- 
sumption C, the transmission losses in any one area are 
not a function of the generation in any other area. 

In order to minimize F subject to the constraint vec- 
tor /, the following set of homogeneous partial differen- 
tial equations are solved: 



SPE ui 



SH 



* Xa + \a 



SPt„ 



10 



15 



20 



25 



30 



35 



SPE^ 



SH 



*X„ + X„ 
X 2 = X^. 

- : K + K 
x a 



&PE„ l 
6PL„ 



X, = 0, or 



SPL U 



SPE ni 



y r = 0, or 



(»•:&-) 



«0, or 



For Company B, 



sh 



SH 



+ X l( - X l =0, or 
SPU 



SH 



-v(»+-^-) . 

( SPU \ 

( I+ splL ) 



&PE cl 



\*PE fl } 

( 6F£ « \ 



hPLu 

K — X a — 0, or 



SH 



PG, 

SH 
PE> 



For the illustrated system, we have 

(a) /„ - PD a + PL a + (PE nl + PE M + PE tt3 ) - PG a 
0 

= PDt + PI* + (PEti + PEm) — PGb = 0 

f e = PD < + PL* + (P^,) - PC, ~ 0 
/, = />£„, +/»£>,«0 

f^PE^ + PEto^O 

f»~PE ta + PE cl <=0 ■ . ■ 

(b) 



Inspection of the equation set shows that the multipli- 
ers associated with the (first three) constraint equa- 
tions/,, above, represent the conventional effective in- 
40 cremental cost of power delivered to their own load 
from their own generation in each of the (three) com- 
panies, while the remaining («) multipliers represent 
the effective incremental cost of power delivered to the 
several tie-lines. 
45 The (three) equations (b) are independent and indi- 
vidually yield the usual intra-company penalty factors 



50 



d lL- / ( x - §W^V 
dGi / I 1 SPG, J 



55 




SH 
SPG b 

SH 

SPG t 



SF b 

SPG t +A * SPG* 



At: this juncture, consideration of an operating re- 
quirement simplifies the problem and facilitates further 
analysis. In order for member companies of a pool to 
benefit mutually from the economies achieved through 
overall, economic dispatch, an equitable method for 
sharing those savings is required. Since the flow of 
power across individual ties cannot readily be con- 
60 trolled, the usual methods for pricing interchanged 
power are based on the average incremental cost of 
power at the boundaries of a given company. Accord- 
ingly, it is reasonable to adopt as an additional criterion 
for economic dispatch the requirement that the aver- 
age incremental cost of power at the boundaries of the 
several member companies be, equal. This criterion is 
in accordance with, but is not required by, assumption 
D. 



65 
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&PE„i &PE„ 3 &PE, a 

X) bPE„ +X « &PE„ + Aa &PE„ 



in area B 



and in area C 



6PE ct 
*• 8P£ f 



* X3 = An , 



10 



We can now calculate such average incremental 
costs of power at the boundaries. A weighted average 
incremental cost, X H , is calculated as shown below, it 
being noted that the sum of the partials for each com- 
pany equals 1. For the 3 company case illustrated, 
in area A 



10 



IS 



20 



since, in the case of company C, power exported out- 
side the pool, or interconnection system, is not to be 
considered. From the above, it is seen that while the in- 
cremental cost of power delivered across the different 
tie lines (X,, X 2 , X 3 ) differs, the weighted average tie 
flow incremental cost for each area is constrained to be 
equal to a common reference lambda, A*. The partials 
of the individual tie flows with respect to net tie flow 
represent participation factors which can be obtained 
as shown below. 

Substituting in these expressions for average incre- 
mental costs at the boundaries the expressions for indi- 
vidual incremental costs at the ties given by equations 
(c), we have for area A 

*• [( ,+ t&r) 6 £t~ + ( 1+ ife-) 



IS 



30 



35 



block diagram form the primary calculation steps in 
solving for the final dispatch signals. Each of the four 
steps illustrated in FIG. 3 is outlined further hereinbe- 
low, and with sufficient detail to allow a programmer 
of ordinary skill in the art to program such calculations 
on a general purpose digital computer having a stored 
program. Steps 71 and 72 are illustrated in detail in 
FIG. 4; steps 73 and 74 involve standard computer op- 
erations well known in the art. 

71. Calculate participation factors. 

For each company (company A is illustrated): 

1. Read the tie flows: PE nl , PE, j2 , PE, i3 

2. Sum the tie flows to form the net tie flow: 

PE„~PE„i + PE M + PE ni 

3. For each such tie, take the difference between the 
current tie flow and the next most recent tie flow, to 
obtain the change in tie flow: 

MC, (') - PE« ) = awe*; 
PE<« (O - PEat {t-l ) = U>E al ; 
PE ai O) ~ PE^ (f-1 ) = 

4. Take the difference between the current net tie 
flow and the next most recent net tie flow to give the 
change in net tie flow: 

5. For each tie, divide the change in tie flow by the 
change in net tie flow, to give the per unit change in tie 
flow defined as CHI(t): 



CHI„(t) m 



&PE„ 



for tie j in company A. 

6. Filter the per unit changes recursively, for each tie, 
by solving the following set of equations, thereby ob- 
taining an estimated per unit change in tie flow defined 
in ECHl(t): 

AU) = 8M(rl)+C«/(()] A(O)=0 
Bit) m 0[B(t-\) + C(t-\) + e*CH!{i-l)] BiO) 
= 0 

C(/) m 0[Cit-\) + 8*CHHt-\)) C(O)«0 
ECHiJLt) m M<r) (J -*)-*<>(!-*)]/* 
ECHI(t) is calculated, from the above cited equations, 



&PE*t 



&PE„ 



Rearranging, and recognizing that the participation 
factors for any one company must add to unity, we 
have 



(■♦a-)-.- 

( ,+ B&-)- 



; K, * CF n 



X„;X 6 *C/>=X, 



for each tie line in each company. The value of 0, as 
used in these equations, is a fraction of about 0.98, 
picked empirically from test runs. It is to be noted that 
50 this is one manner of obtaining the estimated CHI func- 
tion, or ECHI(t), and that alternate means of estimat- 
ing can be used. 
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It is interesting to note that the same set of relation- 
ships are obtained if, instead of imposing constraints /j 
as shown above, on the individual ties, there had been 
imposed one constraint on the sum of all the tie flows, 
which would be a constraint on the sum of the net tie 
flows. Thus, for the three company illustration, the 
fourth constraint equation would be /* = PE„ + PE b + 
PE C = 0. In such case, the one associated multiplier 
(X„) turns out to be identical with the common average 
incremental cost at the boundaries, as defined above. 

Referring now to FIGS. 3 and 4, there are shown in 
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7. Sum the filtered per unit changes in tie flows: 



J=3 

X ECHI4 it) « ECHl nX it) + ECHU it) + ECHU 0) 
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8. Calculate the participation factor for each tie line 
interconnecting the given area (A) by dividing filtered 
per unit changes in the tie flows by the sum of the fil- 
tered per unit changes: 



P.F^ = 



fiPE„j ECHl^it) 
6PE. B ZECHltiit) 
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The above procedure is then repeated for each area control error) signal, (and, if necessary, the system 

(company) within the interconnection system. variables required to calculate the ACE signal), and 

72. Calculate compensation factors for each com- add the total system generation signal and the ACE sig- 
pany, nal to obtain total system load. See the co-pending ap- 

1. Read into the system computer the incremental 5 plication, Ser. No. 163,894, of the same inventor, titled 
transmission losses COMPUTER CONTROLLED COORDINATION OF 

REGULATION AND ECONOMIC DISPATCH IN 

(POWER SYSTEMS, for a discussion of this manner of 

*£ki-S deriving the total system load. 
5P£ " / .10 2. From the total adjusted system incremental cost 

. ■ . , « * curve, determine the system running cost (X R ) corre- 

for the tie lines^ either from on-line load flow or from di tQ nt , oad 

the B-matnx. The procedure for utilizing either on- hne \ Fr * m the adjusted company incremental cost 
flow or the :B-mat™ .is well known ia the art. If on-line cu determine for ^ ch respe ctive company the 
load flow is used, the incremental transmission losses generation to be made available at ^ determined run- 
are pbtamed directly. If the B-matnx constants are ^ CQStj or the ^ d generations: 
used, then the system computer denves the transmis- po o; po^ pg c 

sion losses from the matrix constants. See, for example, 4. Compute the sum of company generation, IPG,, 

O. I. Elgerd: Electric Energy Systems Theory, McGraw and compare with the total system load (PD); If the 

Hill, N.Y., 1971, pp. 294-299, relating to on-line load 2Q sum of company generation is not equal to such total 

flow; N. Cohn: Control of Generation and Power Flow load, adjust the calculated company generations pro- 

on Interconnected Systems, Wiley, N.Y., 1966, pp. portionately. Then, 

71-78, relating to the B-matrix technique. For companies that cannot receive a generation, or 

2. For each tie line within an area, multiply the incre- megawatt signal, 

mental transmission loss for that tie by the participation 25 5. Divide the system running cost (X fl ) by the respec- 

factor for that tie, and take the summation of all ties at tive company compensation factor, to determine the 

the boundary of such company, to yield a weighted av- company running cost (X/). 

erage incremental net tie flow transmission loss for 6. Transmit to each such company its assigned run- 
such company. ning cost. 

3. Calculate the compensation factor for each such 3Q From the above, it is seen that the component curves 
company by adding unity plus the incremental net tie for each individual company are accurate only in the 
flow transmission loss for that company. neighborhood of the current operating point, i.e., the 

73. Calculate the effective interconnection system existing generation and existing actual tie line flows. To 
(pool) cost curve. the extent that the existing generation does not satisfy 

1. Read into the system computer the incremental 35 the economic dispatch criteria, the solution vector (the 
cost curves for the member companies, as transmitted company lambdas) obtained by the algorithm will in- 
from the respective member company computers. It is volve coordinates more or less distant from the current 
to be noted that such curves are transmitted in digital operating points of the companies, thus introducing er- 
form, i.e., a series of data points. For each company, ror. However, assumptions A and 6, together with suf- 
; the data is read in within a range (e.g., 1 0 MW) around 4Q ficiently frequent recourse to the algorithm, will ensure 
.the current operation point. It is assumed that available convergence to an accurate solution vector. Each solu- 
' power is limited to such range, and the company will tion of the dispatch algorithm is based on the values of 

not be asked to deliver power outside of the range. the variables most recently observed or estimated. If 

2. Multiply the company incremental cost curves by the time interval between solutions is small enough, the 
their respective compensation factors, i.e., multiply the 45 change in these variables will be within the limits of ac- 
incremental cost for each data point on the transmitted curacy of the other data used by the solution, 
company curve by the compensation factor for each It will be noted that, as a practical matter, the method 
company. and system of this invention will often be applicable 

3. Determine the minimum and maximum point for only to steam generating units, as hydro units have rela- 
each incremental cost curve, and from these determine 5Q tively invarient running costs. For this reason, step 74 
the most maximum of the maximums and the most min- ( as illustrated in FIG. 3 ) is shown to include reading the 
imum of the minimums, which latter values delimit a pool steam generation. However, it is to be understood 
range of incremental cost, R(X). Each cost curve hav- that the invention is generally applicable to all types of 
ing a maximum or minimum within the limits of R(X) generators without limitation, even though in practice 
is straight line-projected from such maximum or mini- 55 certain operating generators may be excluded from 
mum at a constant MW value. control as described herein. 

4. Subdivide R(X) into an arbitrary number of parts, Referring now to FIG. 5, there is seen an illustration 
e.g., 100, and determine the generation available from of means provided for controlling actual generation at 
each company at each such part, or increment. the company, or sub-area level. Each generator 62 

5. Sum the generation available from all of the com- 6Q (designated as Q lmt where i represents the company 
panies within the system available at each increment, and m represents the number of generators in each 
to yield the composite adjusted curve of total genera- company) has incorporated thereinto a conventional 
tion available at each cost increment. control unit 63 designated as C, m , and adapted to re- 

74. Calculate company power (MW) assignments ceive control signals and to control the power output 
and running costs (X). 65 of the generator 62 in accordance therewith. It is seen 

1 . Read into the system computer generation data . that data representing the power output of each gener- 

from each company, and sum same to obtain a total ator, as, for example, determined by a wattmeter and 

system generation signal. Read the system ACE (area translated into digital form by a conventional sampling 
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and A-D equipment, is transmitted to the respective 
company computer 60, In addition, following the deter- 
mination of the system dispatch signals by the system 
computer 50, and transmission of a corresponding 
company signal back to each company computer 60, 5 
each company computer 60 calculates a megawatt sig- 
nal for each generator 62 within the company, and 
transmits same to the control device 63. As stated here- ' 
inabove, the calculation of the company level dispatch 
signals may be performed in any prior art manner. 10 
Methods of calculating dispatch signals at the company 
level are well known in the art, and consequently are 
not set forth in detail here. 

FIG. 5 also illustrates the use of devices 64 for deter- 
mining the tie flows and generating therefrom digital 15 
signals in proper condition for transmission to the re- 
spective company computers 60. Each device 64 may 
be a conventional wattmeter, along with a sampling 
unit, an analog filtering device and an analog to digital 
converter, for generating periodic representations in 
digital form of the actual tie flow. 

From the above it is seen that there is provided a 
means and method for periodic explicit calculation of 
a dispatch control signal, which calculation is based 25 
upon equalization and optimization of the running cost 
of power delivered to the boundaries of each company 
(or sub-area) within an overall interconnected system. 
The unique algorithm presented hereinabove is not it- 
erative, but enables speedy computer calculation of the 30 
reference system running cost, from which each com- 
pany computer can make accurate and independent re- 
calculations for dispatch of power generation within its 
own company. It is to be noted that, once each com- 
pany is given an assignment (or dispatch signal) opti- 35 
mized in terms of tie flow cost, each company can use 
its own criteria in dispatching generation within its 
boundaries, without upsetting the system optimization 
or the manner of dispatch in the other pool companies. 

I claim: 40 

1. An interconnected power system comprising: 

a. a plurality of subsystems each having a plurality of 
power generation units coordinated by a local gen- 
eration coordination control system for developing 
generating requirements for said units in response 45 
to subsystem generating requirements, said subsys- 
tems being interconnected to transmit power 
among each other; 

b. system digital computer means having a stored 
program for performing, at least at the rate of sig- 50 
nificant change of subsystem generating condi- 
tions, the function of explicitly calculating, for 
each of said subsystems, a subsystem dispatch sig- 
nal as a function of power flow within each such 
subsystem, and 55 

c. transmission means for coupling, at least at said 
rate, status information from said local generation 
control systems to said system computer means, 
and for coupling, at least at said rate, correspond- 
ing subsystem dispatch signals back to control and 
to coordinate the operation of said power genera- 
tion units. 

2. An interconnected power system comprising: 

a. a plurality of subsystems, each including a plurality 65 
of power generating units, said subsystems being 
interconnected to transmit power among one an- 
other; 



b. means for sampling current load and tie line flow, 
for each of said subsystems, at a frequency at least 
as great as the rate of significant change of said cur- 
rent load and tie line flow quantities; 

c. means, responsive to said means for sampling, for 
explicitly developing for each of said subsystems, at 
least at said ^frequency, a subsystem dispatch re- 

.. quirement; 

; d. a plurality of subsystem control means, corre- 
sponding respectively to said.subsystems, for devel- 
oping power generating requirements for the units 
of the associated subsystem; and 
e. means for coupling dispatch requirements from 

. said means for developing, at said frequency, the 
subsystem dispatch requirements to the corre- 
sponding subsystem control means. 

3. An interconnected power system comprising: 

a. a plurality of subsystems, each including a plurality 
of power generating units, said subsystems being 
interconnected to transmit power among one an- 
other; and 

b. system control means including: 

1 . means for sampling prevailing power generating 
conditions at the subsystem level, 

2. principal control means, responsive to the sub- 
system level conditions from said means for sam- 
pling, for explicitly developing, for each of said 
subsystems, a subsystem dispatch requirement at 
a rate as least as frequent as substantial system 
conditions change, 

3. a plurality of subsystem control means, corre- 
sponding respectively to said subsystems, for de- 
veloping power generating requirements for the 
units in the associated subsystem, and 

4. means for coupling dispatch requirements from 
said principal control means to the correspond- 
ing subsystem control means at at least said rate. 

4. The system as described in claim 1, wherein said 
computer means performs said function subject to the 
constraint that the incremental costs of tie flow at the 
boundaries of each such subsystem are equal. 

5. The system as described in claim 4, comprising 
control means in communication with said computer 
means, for controlling power generation units within 
each said subsystem as a function of said subsystem dis- 
patch signals. 

6. A system as described in claim 2 wherein said 
means for explicitly developing comprises computer 
means for developing a subsystem dispatch require- 
ment subject to the constraint that the incremental cost 
of tie flow at the boundaries of each subsystem are 
equal. 

7. A system as described in claim 3 wherein said prin- 
cipal control means comprises system digital computer 
means having a stored program for performing the 
functions of: 

1 . calculating participation factors for each tie line of 
each said subsystem; 

2. calculating as a function of said participation fac- 
tors, compensation factors for each said subsystem; 

3. adjusting the incremental cost curve of each such 
subsystem by its respective compensation factor, 
and summing said adjusted incremental cost curves 
to obtain an effective system incremental cost 
curve; 
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4. determining the system incremental cost from said having N tie lines connected thereto is defined by the 
effective cost curve for the current system load; relation 

and 

5. for at least one of said subsystems, determining a 
generation assignment signal On the basts of the 
system incremental cost and the adjusted cost 
curve of such subsystem, said generation assign- 
ment signal being a form of economic dispatch sig- 
nal suitable for regulating the generation of said 

subsystem. 10 10 - The system as described in claim 9, wherein said 

8. The system as described in claim 7> wherein said participation factors are 
computer means performs, for a second one of said 

subsystem, the function of dividing said system incre- 
mental cost by the calculated compensation factor of 5££n_ — . 
such second subsystem to determine an incremental 15 BPEl * fiA>£l BPE ' 
cost signal for said second subsystem. 

9. The system as described in claim 7 wherein the and add to unity, 
compensation factor (C.F.) for an operating area (i) ***** 
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